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We consider the possibility of solar neutrino decay as a sub-leading effect on their propagation 
between production and detection. Using current oscillation data, we set a new lower bound to 
the 1^2 neutrino lifetime at T 2 / m 2 > 7.2 x 10“^ s.eV“^ at 99%C.L.. Also, we show how seasonal 
variations in the solar neutrino data can give interesting additional information about neutrino 
lifetime. 


INTRODUCTION 


Beyond any reasonable doubt, it is now established 
that neutrinos have non-zero, non-degenerate masses 
and, thus, it would be possible - if not mandatory - for 
them to decay into other particles. 

Before the establishment of the LMA-MSW solu¬ 
tion [H, decay was studied both by itself and in com¬ 
bination with standard flavor oscillations to explain the 
difference between the expected solar neutrino flux from 
nuclear fusion processes in the Sun and the detected flux 
on Earth - the so-called Solar Neutrino Problem (SNP). 

Although it is now ruled out as a leading process Q, 
one can investigate neutrino decay as a sub-leading effect 
in the propagation of solar neutrinos and set limits to 
their lifetime using the most recent experimental data. 

Solar neutrinos are produced in the nuclear fusion pro¬ 
cesses that power the Sun. In such processes, Hydrogen 
nuclei are converted into Helium through several interme¬ 
diate reactions, among which some produce neutrinos in 
very particular spectra - both continuous and monochro¬ 
matic. 

Over the years, several experiments were developed 
for the detection of solar neutrinos at different energy 
ranges. From the pioneer Homestake [H chlorine experi¬ 
ment - which first hinted at the SNP - through the gal¬ 
lium experiments GALLEX Q, SAGE [H and GNO @ 
to the water Gherenkov detectors Kamiokande, Super 
Kamiokande Q and SNO Q. Most recently, the Borex- 
ino 0 experiment also measured the ^Be line. 

The LMA-MSW solution, now considered the best ex¬ 
planation for the SNP, in combination with the measure¬ 
ment of the other oscillation parameters by experiments 
designed for atmospheric, reactor and long-baseline neu¬ 
trinos established the scenario of three massive light neu¬ 
trinos that mix [l|. With such precise measurements of 
the standard oscillation parameters, it is possible to in¬ 
vestigate new phenomena such as the neutrino decay sce¬ 
nario: v' ^ v -\- X. 

For solar neutrinos, the decay of the mass eigenstate 
1^2 into the lighter state vi is disfavored by the data 


and the current bound to z /2 lifetime for invisible non- 
radiative decays [H is r 2 /m 2 > 8.7 x 10“^ s.eV~^ at 
99% G.L.. Similarly, from the combined accelerator 
and atmospheric neutrino data, lifetime is rs/ms > 
2.9 X 10“^^ s.eV“^ at 90% G.L. |iq|. Also, an analy¬ 
sis of long-baseline experiments MINOS and T2K give a 
combined limit of rs/ms > 2.8 x 10“^^ s.eV“^ at 90% 

c.L. [iH. 


In this work, we consider the decay scenario in which 
all the final products are invisible. We combine the avail¬ 
able solar neutrino data with KamLAND [l2| and Daya 
Bay 0 data. For both experiments, the effect of neu¬ 
trino decay is minimum thus allowing us to constrain the 
standard neutrino mixing parameters independently of 
the decay parameter r 2 /m 2 and leading us to obtain a 
robust bound on 1^2 lifetime. Additionally, we show how 
seasonal variations in the solar neutrino data, which are 
enhanced by neutrino decay, can give some interesting 
information about neutrino lifetime. 


FORMALISM 

After production in the solar core, neutrinos propagate 
outwards undergoing flavor oscillation and resonant fla¬ 
vor transition due to the solar matter potential. After 
emerging from the solar matter, they travel across the 
interplanetary medium until they reach the Earth’s sur¬ 
face where they can be detected promptly (during the 
day) or after traversing Earth’s matter (during the night 
- on which they may also be subject to matter effects). 

The transition amplitude for an electron neutrino pro¬ 
duced in the Sun to be detected on Earth as a neutrino 
of flavor a, Vq,, for the standard case of neutrino 

oscillations with MSW effect, can be written as 0: 

A,^=Y^A%AT ( 1 ) 

where is the transition amplitude of an electron neu¬ 
trino produced in the solar core to be in a vi state in 
the solar surface, is the propagation amplitude be¬ 
tween Sun and Earth surfaces, and Af^ is the transition 
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amplitude of a Vi to be in a Vq state upon detection on 
Earth. 

Considering the current limits to their lifetime, neutri¬ 
nos do not decay inside the Sun. It suffices to consider 
their decay on the way to Earth by taking the decay sur¬ 
vival probability for invisible decay, of a neutrino mass- 
eigenstate 2 , with energy after propagating a distance 
L, to be: 


= exp 



with ai = 


( 2 ) 


where rrii is the eigenstate’s mass , is the eigenstate’s 
lifetime and L is the Sun-Earth distance. 

Eor the assumption that only the V 2 mass-eigenstate 
is unstable, the electron neutrino survival probability in¬ 
cluding decay and oscillation is: 


P{Ve ^e) 


4 P?iP?e+Pf2 P2e 


+ 4 , ( 3 ) 


where Sij = sin Oij and Cij = cos Oij and is given in 
Eq. ([2]). In this scenario, one interesting point is that the 
sum over all probabilities is not equal to 1, as explicitly 
we have: 
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Figure 1. Allowed regions for the decay parameter r^lvn^ and 
the mass squared difference Am|i. The hollow curves repre¬ 
sent the solar neutrino analysis only and the filled curves rep¬ 
resent the combined analysis of solar, KamLAND and Daya 
Bay data. The dotted, dashed and continuous line represent 
respectively 90% C.L., 99% C.L. and 99.9% C.L.. 


^ P{V, = C?3 P® (1 - P2^-) . (4) 

Q;=e,/i,,r 

This non-unitary evolution was discussed in Ref. 

Eor this scenario, we should compute P{i'e —^ ^e) and 
^ ^a) that are independent probabilities. 
One important point is that, for appreciable values of 
T 2 /m 2 , the solar neutrino data can be explained by a 
combination of standard three neutrino MSW oscillation 
and decay, which leads to a degenerescence between neu¬ 
trino parameters, specially Am^i and r 2 /m 2 [H . 


ANALYSIS AND RESULTS 

Eor the analysis of z /2 decay over the Earth-Sun dis¬ 
tance and how it affects the expected rate for each solar 
neutrino experiment, we calculate the neutrino survival 
probabilities as shown in Eq. ([3|) and Eq. (|4]), numeri¬ 
cally, under the assumption of adiabatic evolution inside 
the Sun 0- Then, we compute the expected event rate 
for each relevant experiment and compare it to their data. 

We include Homestake total rate Q , GALLEX 
and GNO combined total rate 0 , SAGE total 
rate [HI, SnperKamiokande I full ener^ and zenith 
spectrum [1^, SNO combined analysis [8| and Borex- 
ino 192-day low-energy data 0 . Then, we build a 
function as a function of the relevant parameters 
Xq = xKtan^ 6*12, Am|i,sin^ 0i3,T2lm2). 

We can add complementary information from the reac¬ 
tor experiments KamLAND [HI and Daya Bay [HI and 
their detection of Ue oscillations. One important point 


that led us toward this analysis is the fact that these re¬ 
actor experiments give precise constraints on Am^i and 
sin^^is. KamLAND and Daya Bay have typical base¬ 
lines of LjEj, ^ 10“^^ s.eV“^ and ^ 10“^^ s.eV“^ re¬ 
spectively. Eor the currently allowed values of r 2 /m 2 , 
one has that P^^^ ^ 1, which implies that, in the context 
of these experiments, decay can be neglected and the rel¬ 
evant neutrino probability is the standard three neutrino 
expression P{i'e ^e) = Pee^ with: 

Pee = 1 - sin^ A21 - 5'i3 siu^ ArUgg , (5) 


where Sij = sm20ij, Aij = and drafj = — 

m|, and we define an effective mass square difference 


sin^ ArUgg = 0^2 sin^ A 31 -h 5^2 sin^ A; 


32 - 


This implies that the standard neutrino analysis for 
three neutrinos of KamLAND and Daya Bay experiments 
can also be used for decay scenario. In other words, we 
can identify = xL decay o^r analysis for both 

experiments. The KamLAND experiment have provided 
a function for the standard three neutrino scenario 
used in Ref. [l2| and available in table format as a func¬ 
tion of tan^^i 2 , Am 2 i and sin^^i 3 . Eor the Daya Bay 
experiment, the Xdb function is available in table format 
provided in the supplementary material from Ref. 0 as 
a function of Amgg and sin^ 6*13. 

Then, we write the combined function for solar, 
KamLAND and Daya Bay data as: 

X^ = X©(tan^^'i2, Am2i,sin^6»i3,r2/m2) + 


+ XKL(tan^ 6 »i2, Am^i, sin^ Oiz) + 

+ XDB(^4e>sin^ ^ 13 ), 


( 6 ) 
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Figure 2. Ax^ for 1/2 lifetime T 2 /m 2 . The black (red) curve 
shows only solar (combined) data analysis. 

where was defined before and over which we can 

promptly marginalize the From Eq. (ED, we find 
the allowed regions for independent parameters tan^ ^ 12 , 
sin^ 013, Am^i, and By marginalizing over the 

first two, we get the allowed region for the mass squared 
difference Am^i and the decay parameter T^jvn^ as 
shown in Fig. [U where the hollow (filled) regions show 
the results for the solar neutrino (combined) analysis. 

The degenerescence between Am^i and r 2 /m 2 is evi¬ 
dent in the hollow regions of Fig. [H where higher (lower) 
values of Am^i and lower (higher) values of r 2 /m 2 are al¬ 
lowed. High values of Am^i are ruled out in the standard 
neutrino scenario because it leads to spectral distortions 
that are disfavored by the solar neutrino data. On the 
other hand, it could be turned into a viable solution at 
the cost of having lower values of r 2 /m 2 . The inclusion 
of KamLAND and Daya Bay data break this degeneres¬ 
cence due to their precise independent measurement of 
Am 2 i and sin^ ^13 respectively. We can now precisely 
isolate the contribution of the decay parameter T^jvn^. 
The complete marginalization over the standard param¬ 
eters results in the curve shown in Fig. [2] of Ax^ as a 
function of From it, we can extract a lower limit 

to the V 2 eigenstate lifetime: 

T 2 lm 2 > 7.7 X 10“"^ s.eV“\ at 99% C.L., (7) 

which corresponds to the upper bound to the decay pa¬ 
rameter 0(2 < 8.5 X 10“^^ eV^. 

SEASONAL EFFECT 

One interesting consequence of the decay scenario that 
has never been explored is its effect in the seasonal vari¬ 


Experiment 

Cexp i CTexp 

(Cexp ± ^exp) /Co 

Borexino 

0.0398 ±0.0102 

2.38 ±0.61 

SuperKamiokande-I 

[2^ 0.0252 ± 0.0072 

1.51 ±0.43 

SNO Phase I 

0.0143 ± 0.0086 

0.86 ±0.51 


Table L Experimental best-fit values and errors for Earth’s 
orbital eccentricity e for different solar neutrino experiments. 
We also show the ratio between the fitted values and the 
Earth’s eccentricity eo- 


ation of solar neutrino flux. 

In the absence of decay, the neutrino flux arriving on 
Earth is given by 0® = 0 ®/( 47 rr^), where r = r{t) is the 
time-dependent Earth-Sun distance. The ratio between 
maximum (perihelion) and minimum (aphelion) fluxes is 
i?o = (1 + ^ 0 )^/(1 ~ ^ 0)^5 where cq = 0.0167 is the eccen¬ 
tricity of Earth’s orbit. 

The inclusion of decay modifies the ratio between max¬ 
imum and minimum neutrino fluxes and hence also the 
expected eccentricity e as given by: 


V(rmin) ^ {1 + ef 
V(w) (1 - e)2 


(8) 


where rniax(^min) is the aphelion (perihelion) distance 
and N is the number of events calculated from the con¬ 
volution of the adequate probabilities and cross sections 
for each experiment. 

Erom Eq. ([3|) and Eq. (|4]), we know that 
A^(^min) > A^(^max) d^c to dependence on the 

orbital distance. This implies that R > Rq for any neu¬ 
trino energy and thus, for any neutrino decay scenario, 
an enhancement in the seasonal variation of the solar 
neutrino flux would be expected. 

Thus, the measurement of an eccentricity e > eo is 
a hint in the direction of the neutrino decay scenario. 
In fact, some experiments have measured Earth’s orbital 
eccentricity to be different than the standard value albeit 
still compatible with eo as shown in Table ©• 

Eig. [ 3 ] shows the dependence of the neutrino eccen¬ 
tricity e with the neutrino lifetime T 2 /m 2 as it would 
be measured by different experiments. As it can be 
seen, the lower (higher) the energy of the neutrinos, the 
smaller (bigger) is the enhancement in the seasonal vari¬ 
ation. Since the U 2 content in the neutrino flux leaving 
the Sun depends on the energy, one has that for lower 
(higher) energies, there are less (more) U 2 available to 
decay and hence the seasonal variation will be smaller 
(bigger). Also, due to the decay survival probability, the 
lower (higher) the energy, the bigger (smaller) is the life¬ 
time at which the enhancement is maximum. 

We can include the eccentricity data in the analysis 
with a penalty function added to the x^ for each exper- 
iment: X^easonai = (Sexp - e)V('^exp)^- It results in a 
slightly lower value: 


^■2 /m 2 > 7.2 X 10 ^ s. eV at 99% C.L. (9) 
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Figure 3. Left: Experimental values for e/eo. Black lines are 
the best-fit values and darker (lighter) shades are the la (2a) 
ranges as shown in Table (|T}. Right: Dependence of the orbital 
eccentricity e with the neutrino lifetime T 2 / m 2 as it would be 
measured by different experiments - the ^Be line in Borexino 
(BOR) in red, Super-Kamiokande (SK) in green, and SNO in 
blue. 

due to the fact that the current eccentricity measure¬ 
ments and errors will favor lower, already excluded, life¬ 
times, for which the enhancement in the seasonal varia¬ 
tion (and hence measured eccentricity) is higher. 

CONCLUSION 

We know that neutrinos oscillate with non-zero masses 
and mixing angles. Can neutrinos decay? The answer is 
negative from the combined analysis of data of solar neu¬ 
trino experiments and KamLAND and Daya Bay data. 
From our analysis, we have obtained a new upper bound 
to the z /2 eigenstate lifetime r 2 / m 2 > 7.2 x 10“^ s . eV“^ 
at 99% C.L.. which is almost one order higher than the 
previous bound [H at > 8.7 x 10“^ s . eV“^ at 99% 

C.L.. 

Also, we have shown how decay can enhance the sea¬ 
sonal variation of solar neutrino fluxes and how it affects 
the measurement of Earth’s orbital eccentricity. Cur¬ 
rent data is not good enough to improve the constraints 
to neutrino lifetime. Although future experiments could 
certainly improve on the measurement of solar neutrino 
fluxes and thus better constrain neutrino lifetime, the 
analysis of existing data from later phases of, e.g., Super- 
Kamiokande and SNO for its seasonal variation could, 
in principle, already improve such constraints. We urge 
those experimental collaborations to redo their 

analysis with more of the available data. 

Added Note: During the final stages of this work. 


a paper [ 2 ^ was released which argues for 02 < 9.3 x 
10“^^ eV^ at 2(7 (which corresponds to T 2 / m 2 > 7.1 x 
10“^ s.eV“^). Our result at 2g is 02 < 5.5 x 10“^^ eV^ 
(which corresponds to T 2 j m 2 > 1.2 x 10 “^s.eV“^) 
which is similar to, but a more constrained neutrino life¬ 
time than Ref. [22|. 
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